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ABSTRACT 

We analyse in detail the X-ray data of the microquasar Cygnus X-3 obtained 
during major radio flaring episodes in 2006 with multiple observatories. The analysis 
consists of two parts: probing the fast (~ 1 minute) X-ray spectral evolution with Prin- 
cipal Component Analysis followed by subsequent spectral fits to the time-averaged 
spectra (~ 3 ks) . Based on the analysis we find that the overall X-ray variability during 
major flaring episodes can be attributed to two principal components whose evolution 
based on spectral fits is best reproduced by a hybrid Comptonization component and 
a bremsstrahlung or saturated thermal Comptonization component. The variability of 
the thermal component is found to be linked to the change in the X-ray/radio spectral 
state. In addition, we find that the seed photons for the Comptonization originate in 
two seed photon populations that include the additional thermal emission and emis- 
sion from the accretion disc. The Comptonization of the photons from the thermal 
component dominates, at least during the major radio flare episode in question, and 
the Comptonization of disc photons is intermittent and can be attributed to the phase 
interval 0.2-0.4. The most likely location for Comptonization is in the shocks in the 
jet. 

Key words: Accretion, accretion disks - Binaries: close - X-rays: binaries - X-rays: 
individual: Cygnus X-3 - X-rays: stars 



1 INTRODUCTION 

Cygnus X-3 (Cyg X-3) is a well-known X-ray binary (XRB), 
located at a distance of ~ 8-10 kpc ( Predehl & Schmitt 1995 



Predehl et al.||2000[ [rJickey fc Lockman||1990| ) in the plane 



of the Galaxy. Its discovery dates back to 1966 (Giacconi ct 
aL||1967| but the true nature of the system remains elusive 
despite extensive multiwavelength observations throughout 
the years. The system is thought to harbour a black hole 
due to its spectral resemblance to other black hole XRB sys- 
tems, such as GRS 1915+105 and XTE J1550- 564 (see e.g. 
Szostek et al.|2008| [Hjalmarsdotter et al.|2009[), although a 



neutron star cannot be ruled out (e.g. Vilhu et al. 2009 1 



Distinct 4.8-hour orbital modulation, typical of low-mass 
XRBs, is evident in the X-ray ( Parsignault et al.|1972 I and 
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infrared ( |Mason et al.| [l986) lightcurves, and yet infrared 
spectroscopic observations indicate that the mass-donating 



companion is a Wolf-Rayet (WR) star ( van Keerkwijk et al. 
[1992] ), making it a high-mass XRB. Thus, it is quite possible 
that the accretion disk is enshrouded partially or wholly in 



the strong stellar wind of the WR companion ( Fender et al. 
1999} |Szostek fc Zdziarski|[2008l | Vilhu et al.||2009| ) Unlike 
most other XRBs, Cyg X-3 is relatively bright in the radio 
virtually all of the time. It undergoes giant radio outbursts 
(up to 20 Jy, Waltman et al.|[l995 I with strong evidence of 
relativistic jet-like structures ( Molnar et al.|1 988 , Schalinski 



et al.|1995 Mioduszewski et al.|200l" I . Following the discov- 
ery of transient 7-ray emission from Cyg X-3 (Tavani et al. 
Fermi LAT Collaboration et al.|2009 1 , the 7-ray emis- 



2009 



sion process has been attributed to the interaction of the 



stellar wind with the relativistic jet (Zdziarski et al. 2012 
Dubus et aL]|2010 \ 



Modelling the X-ray data of X-ray binaries and micro- 
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quasars often leads to a problem of degeneracy, i.e. distinct 
multiple models fit the observed data equally well (see e.g. 
|Nowak et al.|[2~011| for a recent take on modelling the X- 
ray spectra of Cyg X-l). Therefore, even if an apparently 
good fit is obtained between the data and the model, it 
does not necessarily imply a match between theory and 
physical reality. In order to make sense out of this degener- 
acy we need to take other data dimensions besides spectral 
into account, namely timing and/or polarisation. While we 
cannot yet access the polarisation dimension, X-ray timing 
data is readily available, and several methods have been 
developed to combine spectral and timing analyses (e.g. 
Vaughan fc Nowak|[l997|) . (jRemillard fc McClintock|[2006 1 



classified the accretion states of X-ray binaries based on con- 
currently examining their X-ray spectral energy distribu- 
tions and power spectra. In this paper, we employ principal 
component analysis (PCA, not to be confused with the Pro- 
portional Counter Array onboard RXTE, always referred to 
here as RX TE/PC A) , one of the standard tools of time series 
analysis (see e.g. Jolliffc 2002 for a review) that has recently 
been introduced in the analysis of the X-ray data of XRBs 
by Malzac et al. ( |2006[ > and references therein. Combined 
with the spectral resolving power of X-ray detectors, X-ray 
lightcurves - instead of representing only one energy band 
- can be transformed into a time series of X-ray spectra. 
Analysing these with PCA yield the principal components 
that are responsible for the spectral variation. This imposes 
a second requirement for the X-ray spectral fits, so that in 
addition to fitting the spectra acceptably, the resulting fits 
also have to satisfy the spectral evolution inferred from the 
PCA. This extra requirement reduces greatly, if not com- 
pletely, the degeneracy of simply using the results from the 
spectral fits in determining the emission components of the 
system. Of course, the PCA does not tell us what the ex- 
act physical components are that correspond to the principal 
components, but these can be inferred by e.g. examining the 
variability spectrum. Ultimately, the physical model can be 
deduced through trial-and-error by fitting the X-ray spectra 
and searching for parameters that correlate with the prin- 
cipal components. iMalzac et al. (20061 showed for Cyg X-l 



that the principal components suggested a variation in pa- 
rameters within one spectral model (eqpair, |Coppi|[l999[ 
specifically the compactness of soft seed photons and the 
compactness of the total power supplied to the plasma). Be- 
low we show that in the case of Cyg X-3 the principal com- 
ponents point to a difference in the normalisation of distinct 
spectral components. 



In this paper we present a detailed analysis of the X- 
ray data from RXTE, INTEGRAL and Swift of Cyg X-3 
taken during the major flaring episodes that occurred in 
2006 May-June. In Section [2 we describe the observations 
used in the analysis, Section 3] introduces the PCA in the 
analysis of multi-energy X-ray lightcurves, and Section [4] 
presents the results from the PCA as well as the results 
from modelling the X-ray spectra. The implications of these 
results are then discussed in Section[5] and we present the 
conclusions of this paper in Section [6j 
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Figure 1. The left-hand side panels from top to bottom show the 
one-day integrated RXTE/ ASM, Swift /BAT and Ryle/ AMI-LA 
15 GHz lightcurves from the May 2006 major radio flare. The 
same is shown for the June 2006 major radio flare in the right- 
hand side panels, with instead the RATAN-600 radio lightcurve 
at 11.2 GHz. The top part of each panel shows the epochs of the 
RXTE (black), Swift (orange) and INTEGRAL (magenta) data 
that are examined in this paper. 



2 OBSERVATIONS 

As mentioned in Section [I] the X-ray observations were 
obtained during the 2006 May-June major radio flaring 
episode, that consisted of two major radio flares with peak 
flux densities of 13.8 Jy (15 GHz) and 11.2 Jy (11.2 GHz) 
at MJD 53865 and MJD 53942, respectively (see Fig. [l). 
For RXTE (both RXTE/PCA and RXTE/HEXTE) and 
SVre/t/XRT data we used the consequent observations taken 
right after a flare and for INTEGRAL (INTEGRAL/ JEM- 
X and INTEGRAL/ISGR1) we used revolutions 437, 438 
and 462 (see Tables [T| [5] and The radio data are from 
( |Koljonen et al.|2010| > 

We reduced all the data as described in the sections be- 
low. In addition to extracting X-ray spectra we constructed 
lightcurves from the data that will be used in the PCA. 
The orbital phase of individual segments of lightcurves were 
determined using a cubic ephemeris ( Singh et al.|2002 l. 
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2.1 Swift 

The five Swift/XRT observations used in this paper are tab- 
ulated in Table [I] Each Swift observation covers approx- 
imately 12 hours and consists of several short pointings of 
~1 ks (which probe different phases of the Cyg X-3 orbit) for 
a total of up to 8 ks exposure times. The data were extracted 
using HEASOFT version 6.12 and reduced with the normal 
procedure using xselect 2.4b. The Swift/XRT lightcurves 
for the PCA were extracted from 62 energy bands spanning 
an energy range from 1.1 keV to 8.0 keV, the first two bins 
having energy bin widths of 0.3 keV, the next three with 
energy bin widths of 0.2 keV, and the rest with energy bin 
widths of 0.1 keV. The background-subtracted 50 s binned 
lightcurve of the energy band 4.2-4.3 keV ranges from 0.6 to 
3.7 counts/s with a mean error of 0.15 counts/s. Likewise, 
the lowest band used (1.1—1.4 keV) has mean counts of 0.5 
counts/s with a mean error 0.09 counts/s, and the highest 
band used (7.9-8.0 keV) has mean counts of 0.12 counts/s 
with a mean error 0.04 counts/s. For the energy spectra we 
use the same energy range as for the lightcurves. In addition, 
we group the data eight bins to one and add 3% systematic 
error for fitting purposes. 



2.2 RXTE 

The 18 RXTE observations used in this paper are tab- 
ulated in Table [2] with four longer observations divided 
into shorter segments (namely 91090-03-01-00, 91090-03-02- 
00, 91090-01-01-000 and 91090-01-01-00) resulting in ap- 
proximately 3 ks long exposure times for RXTE/PCA and 
1 ks long exposure times for RXTE/HEXTE. Each obser- 
vation is individually reduced by the standard method as 
described in the RXTE cook book using HEASOFT 6.12. 
The RXTE/PCA lightcurves were extracted from energy 
bands 4-5, 5-6, 6-7, 7-9, 9-11, 11-13, 13-15 and 15-20 keV, 
whereas RXTE/HEXTE lightcurves were extracted from 
bands 20-30, 30-40, 40-50, 50-70, 70-100, 100-140 and 140- 
200 keV. During these observations only the RXTE/UEXTE 
Cluster B (Cll) was available for lightcurve extraction and 
combined with the rocking of the instrument resulted in 
lightcurves that alternate between 16 s and 48 s intervals. 
Therefore we match the lightcurves from RXTE/PCA to 
those from RXTE/UEXTE to obtain consistent lightcurves 
over the whole energy range of RXTE for the PCA. For 
the RXTE/PCA, in the lowest band used (4-5 keV) the 
background-subtracted count rate ranges from 98.8 counts/s 
to 984.1 counts/s with a mean of 442.4 counts/s and a 
mean error of 5.1 counts/s, and in the highest band used 
(15-20 keV) the background-subtracted count rate ranges 
from 17.4 counts/s to 99.8 counts/s with a mean of 51.5 
counts/s and a mean error of 1.8 counts/s. Likewise, for 
RXTE/UEXTE in the lowest band used (20-30 keV) the 
background-subtracted count rate ranges from 0.9 counts/s 
to 17.5 counts/s with a mean of 8.3 counts/s and a mean 
error of 1.2 counts/s, and in the highest band used (140- 
200 keV) the background-subtracted count rate ranges from 
—3.8 counts/s to 17.3 counts/s with a mean of 0.8 counts/s 
and a mean error of 1.0 counts/s. For the energy spectra we 
use the standard 128-channel RXTE/PCA and 64-channel 
RXTE/UEXTE data with the same energy range as for 
the lightcurves. For fitting purposes we add 0.5% system- 



atic error to the RXTE/PCA data. In addition, we group 
RXTE/UEXTE data to a minimum of six sigma significance. 



2.3 INTEGRAL 

The three INTEGRAL observations are summarised in Ta- 
ble]^ All three observations are approximately 100 ks long, 
comprising six orbital periods with small gaps resulting in 
the exposure times shown in Table [3] The observations 
were analysed using the OSA version 9.0. The lightcurves 
for INTEGRAL/ JEM-X and INTEGRAL/1SGRI were ex- 
tracted from the same energy bands as for RXTE with 
time bins of 100 s. The lowest band (4-5 keV) of 
INTEGRAL/ JEM-X has a background-subtracted count 
rate that ranges from —2.5 counts/s to 29.0 counts/s with 
a mean of 8.4 counts/s and a mean error of 1.9 counts/s. 
Similarly, the background-subtracted count rate for high- 
est band (15-20 keV) ranges from —9.3 counts/s to 18.7 
counts/s with a mean of 1.8 counts/s and a mean error of 
1.8 counts/s. Likewise, the lowest band (20-30 keV) of IN- 
TEGRAL /1SGRI has a background-subtracted count rate 
that ranges from —4.1 counts/s to 40.5 counts/s with a 
mean of 15.3 counts/s and a mean error of 3.5 counts/s. 
The highest band (140-200 keV) of INTEGRAL /ISGRI has 
a background-subtracted count rate that ranges from —6.8 
counts/s to 5.9 counts/s with a mean of 0.1 counts/s and 
a mean error of 1.8 counts/s. For fitting the X-ray spectra 
we use the standard 64-channel INTEGRAL/ JEM-X and 
16-channel INTEGRAL /ISGRI data with a similar energy 
range as for RXTE. We added 3% and 5% systematic er- 
ror to INTEGRAL/ JEM-X and INTEGRAL/ISGRl respec- 
tively. 



3 PRINCIPAL COMPONENT ANALYSIS 

The main use of PCA is to find the smallest number of com- 
ponents that are sufficient in describing the data while losing 
the least amount of information in the process. Basically, 
PCA finds patterns in the data in a way that highlights 
the differences and similarities in the data set. In many- 
dimensional cases, where a graphical representation is not 
convenient, it is a particularly powerful analysis tool. By 
identifying the "new coordinates" of the data set (i.e. the 
principal components) where the data points primarily clus- 
ter and ignoring the small scatter of the points around these 
coordinates, the dimensionality of the data set is greatly re- 
duced, defined only by these new coordinates. In order to 
apply PCA to our X-ray spectra we closely followed the 
work of |Malzac et al.| ( |2006[ ) and the reader is highly recom- 
mended to read that paper and references therein for further 
details. Here, we briefly summarise the main points of the 
analysis. 

We start with a number of spectra p measured at times 
ti, t2, . ■ .t p binned into n bins corresponding to energies 
Ei, E%, . , . E n (see Section 2). The spectra are then viewed 
aspxn matrix with coefficients F(t p ,E n ), i.e. the fluxes 
in each energy band at times p. For the PCA to work, one 
needs to subtract a mean flux F(E n ) from the coefficients for 
each energy channel. One can compute the n x n covariance 
matrix out of the data matrix, which states the variances 
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Table 1. Log of Swift observations 



N° 


ObsID 


Date 


MJD 


N° 


Phase 


XRT 








Interval 


Pointings 


Intervals 


Exp. 






[yy/mm/dd] 


[d] 




[ks] 




1 


00053560010 


06/05/12 


53867.579-53867.990 


7 


0.18-0.24/0.51-0.57/0.85-0.91 


5.64 


2 


00053560011 


06/05/13 


53868.443-53868.992 


8 


0.17-0.25/0.51-0.59/0.85-0.92 


7.37 


3 


00053560012 


06/05/17 


53872.008-53872.949 


10 


0.00-0.07/0.34-0.41/0.67-0.74 


7.86 


4 


00053560013 


06/07/26 


53942.034-53942.975 


12 


0.00-0.02/0.06-0.09/0.34-0.41/0.67-0.74 


3.81 


5 


00053560014 


06/07/27 


53943.027-53943.636 


10 


0.00-0.05/0.34-0.39/0.67-0.72 


5.91 



Table 2. Log of RXTE observations 



N£ 


ObsID 






Date 




MJD 
Interval 






Phase 
Interval 


PCA 
Exp. 


HEXTE 
Exp. 


X-ray/ 
Radio 


Radio 
Flux 










[yy/mm/dd] 




[d] 










[ks] 


[ks] 


State 


[Jy] 


1 


91090-02 


01 


-00 


06/05/12 


53867 


006- 


-53867 


040 





980- 


-1.146 


2.8 


1.0 


FHXR 


7.498 


(Ry) 


2 


91090-02 


01 


-06 


06/05/12 


53867 


725- 


-53867 


777 





580- 


-0.766 


3.2 


1.1 


FIM 


6.348 


(Ry) 


3 


91090-02 


01 


-01 


06/05/13 


53868 


052- 


-53868 


093 





219- 


-0.407 


3.2 


1.1 


FHXR 


6.558 


(Ry) 


1 


91090-02 


01 


-07 


06/05/13 


53868 


181- 


-53868 


271 





875- 


-1.312 


5.1 


1.5 


FHXR 


5.519 


(%) 


5 


91 090-02 


01 


-02 


06/05/13 


53868 


315- 


-53868 


354 


o 


531- 


-0.719 


3.2 


0.9 


FHXR 


5.160 


(Ry) 


6 


91090-02 


01 


-05 


06/05/13 


53868 


836- 


-53868 


894 





154- 


-0.342 


3.2 


1.1 


FHXR 


7.174 


(Ry) 


7 


91090-02 


01 


-03 


06/05/13 


53869 


166- 


-53869 


214 





798- 


-0.982 


3.2 


1.0 


FHXR 


6.678 


(Ry) 


8 


91090-02 


01 


-08 


06/05/14 


53869 


558- 


-53869 


599 





761- 


-0.898 


2.3 


0.8 


FIM 


6.316 


(Ry) 


9 


91090-02 


01 


-01 


06/05/14 


53869 


820- 


-53869 


871 





073- 


-0.261 


3.2 


1.1 


FIM-FHXR 


5.745 


(Ry) 


10 


91090-02 


01 


-10 


06/05/15 


53870 


017- 


-53870 


057 





056- 


-0.245 


3.2 


1.1 


FIM 


5.831 


(Ry) 


11 


91090-02 


01 


-12 


06/05/15 


53870 


214- 


-53870 


254 





041- 


-0.229 


3.2 


0.9 


FIM 


4.829 


(Ry) 


12 


91090-02 


01 


-11 


06/05/15 


53870 


603- 


-53870 


649 





007- 


-0.160 


2.6 


0.9 


FIM 


3.643 


(Ry) 


13 


91090-02 


01 


-09 


06/05/15 


53870 


998- 


-53871 


039 





975- 


-1.164 


3.3 


1.1 


FIM 


1.869 


(Ry) 


14 


91090-03 


01 


-00(1) 


06/05/17 


53872 


051- 


-53872 


084 





243- 


-0.406 


2.8 


0.9 


FIM 


1.353 


(Ry) 


15 


91090-03 


01 


-00(2) 


06/05/17 


53872 


112- 


-53872 


149 





545- 


-0.735 


3.3 


1.0 


FIM 


1.275 


(Ry) 


16 


91090-03 


01 


-00(3) 


06/05/17 


53872 


177- 


-53872 


215 





874- 


-1.063 


3.3 


0.9 


FIM 


1.208 


(Ry) 


17 


91090-03 


02 


-00(1) 


06/05/17 


53872 


767- 


-53872 


804 





825- 


-1.015 


3.3 


1.1 


FIM 


1.120 


(Ry) 


18 


91090-03 


02 


-00(2) 


06/05/17 


53872 


832- 


-53872 


870 





153- 


-0.342 


3.3 


1.1 


FIM 


1.120 


(Ry) 


19 


91090-03 


02 


-00(3) 


06/05/17 


53872 


898- 


-53872 


935 





481- 


-0.670 


3.3 


1.1 


FIM 


1.120 


(Ry) 


20 


91090-01 


01 


-000(1) 


06/07/26 


53942 


459- 


-53942 


499 





827- 


-1.028 


3.5 


1.0 


FHXR 


2.797 


(RA) 


21 


91090-01 


01 


-000(2) 


06/07/26 


53942 


525- 


-53942 


565 





155- 


-0.356 


3.5 


1.0 


FHXR 


2.797 


(RA) 


22 


91090-01 


01 


-000(3) 


06/07/26 


53942 


590- 


-53942 


630 





483- 


-0.684 


3.5 


1.0 


FHXR 


2.797 


(RA) 


23 


91090-01 


01 


-000(4) 


06/07/26 


53942 


656- 


-53942 


696 





811- 


-1.012 


3.5 


1.0 


FHXR 


2.797 


(RA) 


24 


91090-01 


01 


-00(1) 


06/07/26 


53942 


729- 


-53942 


761 





178- 


-0.340 


2.8 


0.9 


FHXR 


2.797 


(RA) 


25 


91090-01 


01 


-00(2) 


06/07/26 


53942 


799- 


-53942 


827 





528- 


-0.668 


2.1 


0.8 


FHXR 


2.797 


(RA) 


26 


91090-01 


01 


-00(3) 


06/07/26 


53942 


868- 


-53942 


892 





872- 


-0.996 


2.1 


0.7 


FHXR 


2.797 


(RA) 


27 


91090-01 


01 


-01 


06/07/26 


53942 


937- 


-53942 


959 





222- 


-0.327 


1.8 


0.6 


FHXR 


2.797 


(RA) 



between each n-dimensions. The PCA proceeds to calcu- 
late the eigenvectors, c„, and eigenvalues of the covariance 
matrix. These eigenvectors form the new coordinates of the 
data and the accompanying eigenvalue states the propor- 
tion of variance of a particular eigenvector, i.e. the highest 
eigenvalue and accompanying eigenvector is the first prin- 
cipal component of the data set etc. One can then form a 
linear decomposition C n k = (c„ic„2 • ■ • c n k) of the data set 
using these eigenvectors, ordered by the proportion of vari- 
ance. The components producing only a small fraction of the 
overall variance can then be dropped, reducing the dimen- 
sionality, i.e. choosing a small number of eigenvectors for the 
linear composition (see more about choosing the number of 
components in Section 3.1l. The data set defined by these 



new coordinates can be obtained by 



Fk P — Ck 



X F n p ^ CkiF%i 



(1) 



where the chosen k coordinates lie in the direction of the 
calculated eigenvectors Ck- Then one can work backwards 
producing the original data set without the less significant 
components, 

k 

Fnp CA = (C nh x F kp ) + F(E n ) = F(E„) + C nl F lv (2) 

i=l 

which most likely are just systematic errors and noise, 
leaving only the intrinsic variability of the data set. Just as 
easily and more importantly, one can concentrate on just 
one principal component, so as to isolate the effect caused 
by this component to the overall variability. In the following 
section we will exploit PCA in different ways using both 
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Table 3. Log of INTEGRAL observations 



N° 


Rev. 


Date 


MJD 


JMX1 


ISGRI 








Interval 


Exp. 


Exp. 






[yy/mm/dd] 


[d] 


[ks] 


[ks] 


1 


437 


06/05/12 


53867.63-53868.90 


75.4 


60.9 


2 


438 


06/05/17 


53872.00-51873.10 


67.4 


55.0 


3 


462 


06/07/26 


53942.25-53943.50 


52.9 


64.2 



spectral and timing information. In the spectral domain this 
includes producing average, minimum and maximum energy 
spectra out of individual principal components (Section 
their fractions of variance as a function energy (Section 



3.21, 



3.31, 



and their correlation with the measured flux. In the timing 
domain we can follow the individual principal components 
by tracking the evolution of Fk p with time t p , for each k = 
1, 2, 3 . . . Similarly, we can test the correlation of Fk p with 
the measured flux with time t p for each energy band E n . 



by dividing the individual variance spectra by the total vari- 
ance spectrum. 



4 RESULTS 

We divide this section into two parts, first concentrating on 
the PCA results from Swift, INTEGRAL and RXTE, and 
then on the fitting of the RXTE spectra. 



3.1 Log-Eigenvalue Diagram (LEV) 

The usual tool for deciding how many principal components 
one should retain is called the scree graph or log-eigenvalue 
(LEV) diagram (see e.g. Jolliffe 2002 for a review). In the 



scree graph, the eigenvalues l n are plotted against n and 
the critical point is where the eigenvalues below it are steep 
and those after it are flat, thus forming an "elbow" in the 
graph. An alternative method - and the one used in this 
paper - is the LEV which plots the logarithm of l n against 
n. If the "noise" components are decaying in a geometric 
progression, the corresponding eigenvalues will appear as a 
straight line in the LEV diagram. Also, instead of plotting 
the eigenvalues l n , we use the proportion of variance, which 
is defined as tk = Ik/ 53 =1 '»> w hh k = 1,2, ... ,n. 



3.2 Extremal energy spectra 

The contribution of an individual principal component to 
the energy spectrum can be determined by looking at the 
Fn P CA from Eq. [2] as a function of energy E„ for each 
k = 1,2,3... The minimum and maximum effect to the 
energy spectra is calculated by integrating F^ p CA over the 
energy range for each p (and k = 1, 2, 3 . . .), and then select- 
ing the minimum and maximum value which then represent 
the minimum and maximum effect to the shape and normal- 
isation of the energy spectrum. 



3.3 Variance spectrum 

The variance spectrum is a graph that shows the measured 
variance as a function of energy E n . It can be plotted for all 
the principal components, i.e. showing the overall variance 
across the energy range of the data, and it is obtained as the 
diagonal components of the covariance matrix C„ n ■ On the 
other hand the variance spectrum can be calculated for the 
individual principal component as the diagonal of the co- 
variance matrix of F£ p CA from Eq. [2] for each k = 1, 2, 3 . . . 
One can then also calculate the contribution of each princi- 
pal component to the total variance as a function of energy 



4.1 PCA results 

The result of the PCA analysis is depicted in Figs. [2j [3] and [4] 
separately for each X-ray observatory. Fig. [2] shows the frac- 
tions of variance attributed to the principal components. 
Fig.[3]shows the variance spectra for the three (RXTE) and 
the two (Swift/ INTEGRAL) most significant principal com- 
ponents, as well as the sum of the remaining components 
labelled as noise. Fig. [4] shows what is the effect of the two 
or three most significant principal components on the Cyg 
X-3 spectra and in what energy bands the variance occurs, as 
well as how it is correlated to the flux in each energy band. 
Overall, the first principal component in all observatories 
contributes the most to the overall variance (86%, 97% and 
91% for Swift/XRT, RXTE and INTEGRAL respectively, 
see Fig. |2J, corresponds to the change in normalisation in 
the X-ray spectra, and is positively correlated with the mea- 
sured flux and thus with the orbital modulation (Fig. [4] bot- 
tom row). Since this component is modulated along with the 
X-ray fluxes, it means that the emission is a local one and, 
thus, not connected to the radio flux, which does not show 
any distinct modulation. In addition to the first principal 
component, there are also other, less significant, principal 
components that contribute in a small way to the overall 
variability. 

In the following we will review the PCA results of each 
observatory independently. 



4.1.1 Swift 

For Swift, the first two principal components add up to 
nearly 90% of the overall variance with most of the vari- 
ance concentrated in the first principal component (86%). 
Based on the LEV diagram (Fig. [5| , there is a sharp elbow 
after the two principal components indicating the start of 
the noise. Noticeable features in the variance spectrum (Fig. 
[3| are dips in energy bins 1.8-1.9 and 2.3-2.4 keV and an 
increase of variability in the iron line region. The dips corre- 
sponds roughly to the location of some of the strongest emis- 
sion lines in the X-ray spectra (H-like silicon at ~ 2.0 keV 
and H-like sulphur at ~ 2.5 keV) and could be interpreted 
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Figure 2. The LEV diagrams of the principal components from all X-ray observatories (left: Swift, center: RXTE, right: INTEGRAL). 
The panels show the proportion of variance attributed to each principal component coloured corresponding to Fig. [3] Also, the individual 
proportion of variance is labelled for the first two or three principal components that are presented in Figs. [3] and |4] 




Figure 3. The variance spectra of all observations (black solid line with data points) from all observatories (left: Swift, center: RXTE, 
right: INTEGRAL). The lower curves in the figure show the contributions of the principal components «i (red, solid), »2 (blue, dash), 
03 (green, dot) and the remainder of the components totalling to noise and systematic errors (orange, dot-dash). 



as a reduction in variability indicating line production fur- 
ther out from the compact object in the photoionised stellar 
wind. However, the dips could also be due to a low count 
rate and wide energy bins. The decrease in variability in the 
second principal component around ~5 keV as visible in the 
variance spectrum is most likely of instrumental origin since 
there are no known emission lines in that region. To confirm 
this, we performed a similar PCA analysis on a set of simu- 
lated lightcurves based on the best-fitting model as we did 



on the original data (see Section 4.3 1. Both principal compo- 



nents show a change in variability in the iron line region, so 
that a decrease in variability corresponds to times of maxi- 
mum spectra and the increase to times of minimum spectra 
(Fig. |4j upper and upper middle rows). Interestingly, the 
first principal component show just a peak around 6.5 keV 
corresponding to mostly cold iron line, while a sharper peak 
at 6.7 keV corresponding to ionised helium-like iron line in 
addition to 6.5 keV peak is present in the second principal 
component. This implies that there are two regions where 
the iron line forms, where the other (corresponding to the 
second principal component) is more ionised and thus hot- 
ter. Both principal components are correlated with the X-ray 
fluxes (Fig. |4j bottom row), although the second principal 
component slightly less but the correlation is increasing to- 
wards the iron line region. Compared to the principal com- 
ponents of RXTE and INTEGRAL, the variance spectrum 
of Swift appears to not to include a third principal com- 



ponent. The analysis of the simulated data in Section |4.3| 
indicates that the spectral components that were found to 
correspond to the second and third principal components 
seen in the RXTE and INTEGRAL data are clumped to- 
gether in Swift data, most likely an effect due to the small 
size of the dataset. 



4.1.2 RXTE 

For RXTE, almost all the variance (97%) is due to the 
first principal component. Compared to the principal com- 
ponents from the Swift analysis, it is more difficult to es- 
tablish where the noise components commence based on the 
LEV diagram (Fig. |5|. However, some levelling off of the 
components is visible after the third principal component, 
and the fourth principal component dips below 0.1%, a very 
small proportion of the variance and thus most likely con- 
tributing to the noise. The amount of RXTE data compared 
to the other observatories is manifold, and thus we can ex- 
pect the noise level to be lower. We opted to include the 
third principal component, albeit contributing only 0.3% to 
the overall variance, since its shape in the variance spectrum 
(Fig. |3| resembles a reflection spectrum. 

The second principal component contributes 2.3% of the 
overall variance. The variance spectrum of this component 
exhibits a cut-off at ~10 keV. Despite not contributing much 
to the variance overall - 20% at most around 20 keV (Fig. 
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Figure 4. The first two or three principal components {a\, 02 and az) for data from all X-ray observatories. The upper panels show the 
average (dash), minimum and maximum (thick lines) spectra of the principal components. The upper middle panels show the ratio of 
minimum and maximum spectra to the average one. The lower middle panels shows the fraction of variance that the principal components 
accounts for as a function of energy and the bottom panels show the correlation coefficients of the principal components with flux in 
each energy band. 



[4] lower middle row) - the fact that it affects the energy 
spectra most in the region 10-40 keV (Fig. [4] upper and 
upper middle row) explains why it is required to fit the X- 
ray spectra (see Section 4.2 1. Interestingly, this component 
is mildly anti-correlated with the X-ray fluxes. Fig. [4] shows 
that the third component impacts the energy spectra and 
the variance minimally, and it is not correlated with the 
X-ray fluxes. 

In addition to the analysis performed above, also 
"colour" (i.e principal component-component) and phase- 
folded principal component diagrams were constructed using 
the RXTE data. These are plotted in Figs. [5] and [6] respec- 
tively. Fig. [5] shows the first principal component as a func- 
tion of the second principal component calculated from all 
the RXTE pointings. We can see that in most observations 
the increase of both principal components happens more or 
less in a linear fashion, thus signalling a relationship between 
the two. Four observations (pointings 7, 8, 10 and 17) have 
been fitted with a linear regression and highlighted to show 
changes in the diagram. Also marked in Fig. [5] are the X- 
ray/radio states during the observations and a change from 
theFHXR state to the FIM state^can be seen approximately 
at the mean value (0.0) of the second principal component. 
As the radio flare fades, the slope of the fitted line decreases 
and the amount of variability in the second component in- 



1 The state nomenclature follows the classification of l |Koljoneii] 
|et al.p blO:) where the FHXR and the FIM state correspond to 
intermediate-hard and intermediate-soft X-ray spectrum with ra- 
dio flaring, respectively 



creases. Fig. [6] shows the principal components phase-folded 
to the orbit of Cyg X-3. 



4.1.3 INTEGRAL 

Similarly to Swift and RXTE, most of the variance (91%) is 
caused by the first principal component. However, the LEV 
diagram (Fig. [2] right panel) indicates that the components 
level off after the first component, and thus the INTEGRAL 
data does not probe well the minor components. Despite 
this, we include the second component (contributing 3.5% 
of the overall variance) in the analysis as it resembles the 
second principal component from RXTE, due to its similar- 
ity in shape with the variance spectrum (Fig.[3|, as well as in 
the variance fraction (Fig. [4] lower middle row), the energy 
spectra (Fig. |4j top and upper middle row), and exhibiting 
the same correlation behaviour with the X-ray flux (Fig. [4] 
bottom row). 



4.2 Model of X-ray spectra 

As mentioned in Section 1, with the aid of PC A the models 
of the X-ray spectra have to fulfil two requirements: fit the 
X-ray spectra with good statistics, as well as producing a 
comparable variation in the fit parameters (whatever they 
may be) matching the evolution of the principal components. 
In the following, we will go through the individual compo- 
nents of our best-fitting model for the X-ray spectra and 
explain our choices in detail. Since the bulk of the spectra 
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Figure 6. The three first principal components phase-folded 
through all the RXTE data. 



Figure 5. The first principal component (»i) plotted as a func- 
tion of the second principal component (c*2) as calculated from 
the RXTE data. Black and grey points refer to FHXR and FIM X- 
ray/radio state respectively as marked in Table [2] Orange points 
refer to pointings 7, 8, 10 and 17 from left to right with a simple 
linear regression fitted for each pointing. 



were RXTE spectra, we use them as a basis for compari- 
son with the PCA results. Simultaneous fits with the Swift 
and INTEGRAL spectra were used to narrow down param- 
eter ranges. All the data fitting was performed using ISIS 
(Interactive Spectral Interpretation System; Houck|2 002). 



4-2.1 Absorption 

One of the crucial points to understanding the underlying 
physical properties of Cyg X-3 is its complex absorption 
profile. There is increasing evidence that the absorption is 
divided into multiple parts. One of the components is the 
absorption of the interstellar medium in the line-of-sight, 
which is estimated to be 1.5 xlO 22 cm -2 (Dickey & Lock- 



1990 1, a result obtained using the 21 cm wavelength ab- 



sorption. However, when studying the X-ray scattering halo 
( jPredehl fc Schmittj 1995) the absorption is found to be 3.0 
x 10 cm - , clearly in disagreement with the former result. 



In addition, Szostek & Zdziarski ( 2008 ) studied the effect of 



the stellar wind on the emanating soft X-ray radiation from 
the compact object and found that the X-ray spectrum is 
well fitted when including cold clumps in the wind resulting 
in a partial covering factor of 1%. Therefore, it would be 
appropriate to use at least three different absorption compo- 
nents in the fits. However, as is the trend in modelling X-ray 
spectra, the number of free parameters in the models vastly 
increase when taking into account more physically realistic 
scenarios. A balance has to be found between a physically 
realistic model and a phenomenological model, but usually 
the former has so many variables that it will fit anything. 



Whole model _ 

BELM 

Brcmsstrahlung : 

Linemodcl 




100 1 
Energy [keV] 

Figure 7. An example of the FHXR (left panel) and FIM 
(right panel) energy spectrum with the best-fitting model (black, 
solid line) and the individual components (Comptonization: blue, 
dashed line; bremsstrahlung: green, dotted line; iron line: orange, 
dot-dashed line) overlaid and labelled. 



Therefore, one has to strive for physically realistic models 
but endeavour to maintain the degrees of freedom as low 
as possible. To achieve this, we model the X-ray absorption 
with the simplest absorption profile, which consists of a uni- 
form component (phabs) and a partially covered component 
(pcfabs). This model has successfully been applied to Cyg 



X-3 spectra previously in Koljonen et al. 
all ( [2003]), |Szostek fe Zdziarskij ( |2004[ ) and 



(2010), Vilhu et 



Hj almarsdotter 



et al. (20081. We froze the parameter values according to 
fits performed with simultaneous RXTE, INTEGRAL and 
Swift spectra and let only the covering fraction vary when 
fitting only RXTE spectra. 
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4-2.2 Line emission and absorption 

Cyg X-3 harbours a large number of emission lines that are 
especially prominent in the high-resolution energy spectra, 
e.g. Paerels et al. (2000). However, due to the low sensi- 
tivities of the RXTE/PCA below ~ 3.5 keV, only the iron 
line complex is important for the RXTE spectra. |Paerels| 
et al. (20001 showed that the iron line complex consists of 



helium-like and hydrogen-like iron ions (XXV/XXVI) at 6.7 
keV and 6.9 keV, respectively, and cold iron Kq line at 6.4 
keV. However, these lines blend into a single broad iron line 
feature in the RXTE/PCA, and therefore we include only 
one gaussian line component in the model to represent this 
feature as a whole. More importantly, though, these ions 
also absorb hard X-ray photons producing helium-like and 
hydrogen-like iron edges in the energy spectra at 8.83 and 
9.28 keV, which are fixed values in our model only their nor- 
malisation is allowed to vary in the fitting process. Due to 
the inclusion of Swift spectra in simultaneous data fits we fit 
a line model by Savolainen et al. (in prep.) to fit the data. 
These lines most likely form in the stellar wind plasma sur- 
rounding the WR companion that are excited by the ener- 
getic radiation field from the vicinity of the compact object. 



However, the PCA results from Swift (Section 4.1.11 indi- 



cates that the iron lines could form in two locations where 
the other one contains more ionised iron. 



4-2.3 Continuum components 

Based on the PCA results we know that the X-ray vari- 
ability is driven mainly by two components. From previ- 
ous knowledge we know that the X-ray spectra of Cyg X- 



3 are well fitted by hybrid Comptonization (e.g. Koljonen 



et al. (20101; Hjalmarsdotter et al. ( 2008 1) and therefore 



we assume that the first principal component is also cou- 
pled to this model. For the Comptonization model we use 
a relatively new model, belm ( jBelmont et al.||2008[ ), that 
calculates self-consistently the equilibrium electron distri- 
bution taking into an account Compton scattering, syn- 
chrotron emission/absorption, pair production/annihilation 
and Coulomb collisions, assuming a spherical, magnetised 
and fully ionised proton-electron/positron plasma with ra- 
dius R. Akin to the more widely used Comptonization model 



eqpair ( Coppi 1999 1, belm parametrizes energy injection as 
compactness 



I 



(It 



a R' 



(3) 



where L is the power supplied to the plasma, or is the 
Thomson cross section and m e is the electron rest mass. The 
power can be supplied in three forms: (1) non-thermal elec- 
tron acceleration (l n th), with the injected electron distribu- 
tion modelled as a power law (index Pi„j, i.e. n e (7) oc -y 
from energies jmin to jmax, (2) thermal heating of electrons 
(hh) e.g. in Coulomb collisions of hot ions, and (3) external 
radiation (l s ) from a geometrically thin accretion disk with 
a pure black body spectrum of temperature kTtb- In a sim- 
ilar fashion to the electron-photon coupling, the magnetic 
energy density (B 2 /8ty) can be parametrized by 



Ib = 



^R B 



8tT 



(4) 



for tangled fields. The model spectra are tabulated for 
fitting in ISIS with fixed parameters to reduce the overall 
degeneracy of the parameters as well as computation time. 
The total compactness of the injected power that results in 
radiation (I = l nt h + hh + h) can be estimated for Cyg X-3 
assuming an average flux of Fo = 1.3 x 10 -8 erg s _1 cm~ 2 , a 
distance of 9 kpc, a black hole mass of 30 Mq, and a plasma 
size of 30 gravitational radii as: 



-(h) ( 



d 
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(5) 



We consider only purely non-thermal electron injection 
(i.e. l„th 7^ 0, Ith = 0) with electron distribution parame- 



ters fixed at r,. 



1.3 and 7„ 



1000. We 



can therefore vary l a and derive the non-thermal compact- 
ness using I = l n th + Is- Since the electrons are allowed to 
thermalise through Coulomb collisions, this produces a hy- 
brid thermal/non-thermal distribution in equilibrium. The 
hybrid model produced good results previously in |Hjalmars-| 



dotter et al. ( 2009 1 and references therein using eqpair and 



Koljonen et al. (20101 using COMPPS 



However, we find that not all the X-ray spectra studied 
are well fitted by using Comptonization alone, not even by 
allowing some of the fixed parameters mentioned above to 
vary as well. Also, as the PCA results show, a second prin- 
cipal component is required to explain all the X-ray vari- 
ability and therefore we include a second continuum com- 
ponent in the model. As its overall effect in the PCA and 
X-ray spectra is smaller than the first principal component, 
i.e. Comptonization, degeneracies arise, and multiple com- 
ponents will fit the spectra. We found good fits when using 



reflection (reflect, Magdziarz & Zdziarski 19951, multi- 



colour disc blackbody (diskbb, Mitsuda et al. 1984 1, thermal 
bremsstrahlung (BREMSS, |Kellogg et al.j!975 l, or other ther- 



mal Comptonization models (compst, Sunyaev & Titarchuk 
[1980} comptt, |Titarchuk|199"4"l THCOMpT jZycki et al.|19"99f 7 
in addition to BELM. However, only a few of the above models 
produced the variability of the principal components, thus 
reducing the number of suitable models. Table [4] shows the 
chi-squared values of all the models together with correla- 
tions against the principal components and normalisation 
of the continuum components fitted to RXTE pointings 1— 
13. The correlations were calculated using robust methods 



(e.g. Rousseeuw & Leroy 1987) that offer good outlier de- 



tection, so that outliers in the data do not detract from the 
"real" correlation. Robust methods supply ways of detect- 
ing, weighing down (or rejecting altogether), and marking 
outliers that largely removes the need for manual screening. 
In addition to robust correlation we use robust regression, 
its coefficient of determination (R 2 , essentially telling how 
close the data points cluster around the regression line) and 
the F-test probability of a null hypothesis in assistance to 
study the strength and significance of the correlation. The 
robust statistics are implemented using the ROBUST library 
( Wang et al.|2012[ ). 

Based on Table|3] we see that a lot of different setups will 
provide good reduced chi-squared values. However, when 
taking into account the correlations between the principal 
components and spectral model normalisations, four models 
(diskbb, bremss, compst, and comptt) stand out in Table 
[4] All of these have optically thin (r = 0.5) Comptonization 
and some thermal (albeit hot) component. These four mod- 
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els were then fitted to the whole dataset (bottom part of 
Table Hi. Based on these fits, the COMPST model performs 
the weakest and thus if the second principal component is 
due to thermal Comptonization it probably does not arise 
from cold seed photons. The diskbb model has slightly bet- 
ter correlation with the second principal component than 
COMPSTmodel. However, the disk temperature is very high 
(~ 5 keV), thus the radiation most likely is not coming from 
the disc. The possibility of a heated disc is discussed in Sec- 
tion |5.5| Basically, the bremss and COMPTT models produce 
equally good fits and correlations. For the COMPTT model we 
set the seed photons to the same temperature as the BELM 
seed photons and the electron temperature was frozen to 
the values of the diskbb model. With these restrictions the 
optical depth for the thermal Comptonization is very high 
(r ~ 5 — - 10), which, with an electron temperature of 5 keV, 
results in a Compton-y of ~4 thus approaching saturated 
Comptonization. To summarise, the best-fitting model in 
spectral and variability terms has optically thin, rather ther- 
mal, Comptonization dominating the variability throughout 
the X-ray regime and a thermal, rather hot (~5 keV) plasma 
component producing variability in the ~ 10-20 keV regime 
and radiating through Comptonization or bremsstrahlung. 
For the rest of the paper we have modelled the second spec- 
tral component as a thermal bremsstrahlung component (see 
discussion about the origin of the thermal component in Sec- 
tion 5.5 1. 

We would also like to note that a marginal effect (~1%) 
on the variance spectrum is caused by the third principal 
component which, based on the shape in the variance spec- 
trum, is most likely the reflection component from the ac- 
cretion disc. As the reflection bump would cause only a very 
minor adjustment to the X-ray spectra above 10 keV we 
chose not to include it in the model, and we consider the 
iron Ka line from reflection to be blended in the gaussian 
line of the model. However, as discussed in the next section, 
the centroid of the gaussian line follows the third principal 
component so that when the centroid approaches 6.4 keV 
(i.e. the Ka line) the proportion of variability of the third 
principal component increases and vice versa. This evolu- 
tion also gives credibility to the interpretation that the third 
principal component is caused by reflection. 

The best-fitting models for all RXTE spectra are tab- 
ulated in Table |A| The best-fitting models for simultaneous 
RXTE, INTEGRAL and Swift spectra are tabulated in Ta- 
ble |A2| with two examples from the FHXR and FIM states 
plotted in Fig.[7|together with individual model components. 
The implications of these results are discussed in Section [5] 

4.3 Simulating the effect of parameter variability 
on the variance spectra 

To double-check that the above best-fitting model indeed re- 
produces the PC A results presented in Section [4. 1[ we have 
simulated lightcurve-spectra for Swift, RXTE and INTE- 
GRAL data to see how the variability of the spectral pa- 
rameters affects the variance spectrum. This simulation also 
serves as a method to check the effects of Swift/XBT instru- 
ment response as mentioned in Section |4.1.1| The method 
for creating the lightcurve-spectra is as follows: the auxiliary 
response files and redistribution matrix files of a pointing 
were fed into ISIS together with the best-fitting model. To 



reproduce the variability reminiscent of that seen in the orig- 
inal lightcurves, the normalization values of belm, bremss 
and emission lines were multiplied by the values of the cor- 
responding principal components. This step also allows to 
switch between different setups, e.g. first allowing only the 
normalization of one component to vary while keeping the 
others at zero and after that adding more components. Then 
the spectrum is faked according to the same exposure times 
that were used in the actual lightcurves (see Section Af- 
ter faking the fluxes from the same energy bands that were 
used in the original analysis were retrieved. This procedure 
is then looped through to obtain the same length for the 
simulated lightcurves as in the original data (with varying 
random seed for the fake it function). The results are de- 
picted in Figs. |8j [9] and [10] In these fi gures the first panel 
from the left shows the variance spectrum (bottom) and 
LEV diagram (top) when only belm was varying, the sec- 
ond panel shows when belm + bremss were varying, the 
third panel shows when belm + BREMSS + emission lines 
were varying, and the last panel shows the original variance 
spectrum and LEV diagram. 

The simulated data of Swift and INTEGRAL show that 
the second component is a mixture of bremsstrahlung and 
line emission, which could suggest a reflection model. How- 
ever, the simulated (and original) RXTE PCA data (which 
is the best data set) show that this thermal (bremsstrahlung) 
component is distinct from the reflection component. Nei- 
ther do the model fits support the reflection component. In 
addition, the whole X-ray spectrum is rather thermal dur- 
ing the major flare (the hard tail rises later on), so there are 
not too many hard photons available for producing a strong 
reflection component. 



5 DISCUSSION 

Once successful fits to the spectra have been generated, we 
can examine in more detail what the results imply in terms 
of parameter evolution as the major flare decays, as well as 
how the parameters vary along the binary orbit. 



5.1 Parameter correlations 

As required, the normalisations of the different spectral com- 
ponents follow the principal components (see Figs 
1 12[ ) , i.e. the first principal component ol\ tracks the nor- 
malization of the Comptonization model (belm) and the 
second principal component 0:2 the normalization of the 
bremsstrahlung model (bremss). The third principal com- 
ponent 03 most likely corresponds to a reflection model (not 
included in the fits) as was discussed in Section [4] Other in- 
teresting possibly correlated parameters include the Comp- 
tonized soft seed temperature (kTbb,com P ) that appears to 
track the bremsstrahlung normalization, which could indi- 
cate that most of the seed photons for Comptonization are 
produced by bremsstrahlung emission. This interpretation 
is reinforced by the linear relationship seen between the 
first and the second principal components in Fig. [5] and 
also by the convergence of the bremsstrahlung temperature 
(kTbb,brems) and kTbb,com P in Fig. 12 However, there are lo- 
cal minima in the kTbb,com P that are marked with dashed 



and 



vertical lines in Fig. [TT| At the same time local minima are 
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Table 4. RXTE spectral model comparison and corresponding statistical parameters. The model type is of the following (except for 
REFLECT which is a convolution model): PHABS X PCFABS X edge(1) x edge(2) x (egauss + BELM + X), where X is an additional 
component marked in the first column. The second column shows the mean reduced chi-square value and its standard deviation of 
the models fitted to the observations marked in the last column. The third column shows the robust correlation, the coefficient of 
determination (R 2 ) and the F-test probability of a null hypothesis for the normalisation parameter of BELM with the first principal 
component. The fourth column is similar to the third except that the correlation statistics are taken from the normalisation parameter 
of the second spectral component in question (first column) with the second principal component. 







Nbelm vs. c«i 




Nx vs. ct2 






Model X 




Corr. / R 2 / Pr(F) 


Corr. / R 2 / Pr(F) 


Obs. 


DISKBB(f ) a 


1.01/0.24 


0.67/0.74/0.01 




0.55/0.07/0.46 


1- 


-13 


DISKBB(2) b 


0.94/0.16 


0.88/0.73/1x10- 


-5 


0.91/0. 85/6X10" 3 


1 


-13 


THCOMP(f) c 


0.99/0.25 


0.87/0.57/2x10- 


-3 


-0.30/0.36/0.06 


1 


-13 


THCOMP(2) d 


0.99/0.25 


0.95/0.68/8x10- 


-5 


-0.36/0.33/0.09 


1 


-13 


THCOMP(3) c 


0.99/0.18 


0.58/0.33/0.03 




-0.17/0.01/0.69 


1 


-13 


BREMSS 


0.92/0.18 


0.86/0.68/2x10- 


-4 


0.95/0. 75/3X10" 6 


1 


-13 


COMPST f 


0.93/0.17 


0.81/0.59/3x10" 


-4 


0.88/0.59/4X10" 4 


1 


-13 


COMPTT s 


0.94/0.17 


0.93/0.74/2x10- 


-5 


0.88/0.75/10~ 6 


1 


-13 


REFLECT 11 


0.90/0.18 


0.86/0.70/3x10- 


-5 


0.38/0.13/0.24 


1 


-13 


diskbb(2) 


0.99/0.17 


0.94/0.87/2x10- 


-16 


0.64/0. 63/3X10- 3 


1- 


-27 


BREMSS 


1.02/0.20 


0.97/0.89/7x10- 


-15 


0.83/0. 70/7x10"* 


1 


27 


COMPST 


0.99/0.17 


0.95/0.83/5x10- 


-14 


0.52/0.39/0.01 


1 


27 


COMPTT 


1.00/0.17 


0.94/0.87/3x10" 


-15 


0.76/0.68/5X10" 8 


1 


27 



a Colour temperature is set equal to the temperature of the seed photons for BELM. 
b Colour temperature thawed. 

c Seed photon temperature is set equal to the temperature of the seed photons for BELM and the electron Y is frozen at 2. 
d The same as above, but with seed temperature frozen to 50 eV. 
° The same as above, but with electron Y thawed. 

f Electron temperature is equal to colour temperature in model diskbb(2). 

Seed photon temperature is set equal to the temperature of the seed photons for BELM, electron temperature is set equal to colour 
temperature in model diskbb(2), and spherical approximation scheme. 

h Iron abundances, metal abundances and cosine of inclination set to 1.48, 1.58 and 0.5 respectively. 




Figure 8. Simulated variance spectra (bottom panels) and LEV diagrams (top panels) for Swift. See text for details. The first panels 
from the left show the LEV diagram and the variance spectrum when only the BELM normalization was varying. The second panel from 
the left shows the LEV diagram and the variance spectrum when BELM and BREMSS normalisations were varying. The third panel from 
the left shows the LEV diagram and the variance spectrum when all the component normalisations were allowed to vary, and the last 
panel shows the actual LEV diagram and the variance spectrum calculated from the data. 
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Figure 9. The same as in Fig. [8] but for INTEGRAL data. 
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Figure 10. The same as in Fig. [8] but for RXTE data. 



found (in most cases) for «2 , the bremsstrahlung normaliza- 
tion, and the gaussian line centroid, while the external seed 
photon compactness and 03 (reflection) reaches local max- 
ima. These times correspond to phase interval 0.2-0.4 (when 
the increased variability and low-frequency QPOs were iden- 
tified by |Koljonen et al.|2011[ ) and one can also see that the 
gaussian centroid approaches 6.4 keV, indicating increas- 
ing cold iron Ka emission. Thus, it is reasonable to assume 
that this is when the impact of the accretion disc is more 
prominent. Based on these results we can reinforce the con- 



clusion in Zdziarski et al. (2009) that there is a thermal 



plasma component in the system in addition to the accre- 
tion disk/corona combination. It is also evident that the 
phase interval 0.2-0.4 is somehow distinct, and it is in this 
phase that we expect to probe deeper into the system. 



5.2 Orbital variations 

Since Cyg X-3 has a very short orbital period, it is possible to 
investigate possible changes in the X-ray spectra as a func- 
tion of phase and relate them to geometrical properties of the 
source. In order to look for phase-dependent changes in the 
emission components acquired from the RXTE spectral fits, 
we have stacked the values of the fit parameters from each 



observation into so-called phase-circle diagrams (Fig. 13 1. 
Each of these diagrams displays three fit parameters plotted 
in concentric rings, where the angle shows the orbital phase 
of the system and the width of the curve shows the mean 
parameter value from the stacked observations with the min- 
imum value on the inner circle and the maximum value on 
the outer circle of the corresponding ring. The phase-circle 
diagram allows a quick assessment of how the fit parame- 
ters vary along the orbit and how they relate to each other. 
Based on Fig. |13| we see the following: 
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Figure 12. The scatter plot matrix of the selected fit parameters 
and first two principal components with robust correlations (e.g. 
|Rousseeuw &; Leroy|1987| l drawn and written in appropriate grid 
cells. The black lines show the best linear robust regression and 
the grey ribbons show the standard 95 % significance limits. The 
red numbers correspond to the R-squared value (i.e. the goodness 
of the robust fit or the correlation between the predicted values 
with the actual values.) over 0.6. Thus, coloured numbers close 
to ±1 imply correlation. 



Figure 11. Lightcurves of RXTEj ASM, Swift/BAT, Ryle, se- 
lected fit parameters and three principal components during the 
May outburst. The vertical dotted lines mark observations (point- 
ings 3, 6, 9, 14 and 18) where the seed temperature for Comp- 
tonization expresses a local minimum. These times correspond 
somewhat to local minima in the second principal component, 
bremsstrahlung normalization, gaussian line centroid and to max- 
ima of soft photon compactness for Comptonization, and are at- 
tributed to phase interval 0.2-0.4. See text for implications. 



(i) The covering fraction of the partial absorber has 
a minimum at phase interval 0.75-0.00. During this phase 
interval the iron line centroid and the equivalent width ap- 
proach maximum values and the line width approaches min- 
imum value. Thus, during this phase we see a region which 
is not as strongly absorbed and produces a strong iron line 
that is highly ionised but narrow. A possible interpretation 
could be that the compact object has swept up the wind and 
hence is at least temporally in a lower density region; 

(ii) The normalization of Comptonized emission and the 
external photon field compactness follow the overall phase 
variation of the bolometric X-ray flux. This is in line with the 
notion that almost all the X-ray flux observed comes from 
the Comptonized emission of external photon field that in 
turn photoionises a region producing an iron line. The op- 
tical depth of the H-like iron edge and the equivalent width 
of the iron line show anti-correlated behaviour to the X- 



ray flux so that the maximum occurs around phase 0.0 and 
minimum at phase 0.5. This is most likely due to increasing 
stellar wind depth along the orbit; 

(iii) During phase interval 0.3-0.4 the normalization of 
the bremsstrahlung emission and the temperature of the 
soft seed photons undergoing Comptonization is at mini- 
mum. This might be an effect of the dearth of available data. 
However, the phase interval 0.2-0.4 also corresponds to the 
epochs marked with dotted lines in Fig. |11| where the cen- 
troid of the gaussian line complex dips closer to 6.4 keV. It 
is reasonable to assume that during this phase the geometry 
of the system changes, e.g. the disc photons are able to see 
the Comptonizing component, which is also supported by 
the increase in variability and the presence of low-frequency 
QPOs during this phase ( Koljonen et al.|2011 l; 

(iv) The normalization of bremsstrahlung emission, 
which is proportional to the electron and ion densities and 
the volume of the emission region, has two maxima at phases 
0.0-0.25 and 0.5-0.75. This also explains why the second 
principal component as seen in RXTE and INTEGRAL data 
is anticorrelated with the X-ray fluxes Fig. [4] Therefore, ei- 
ther the electron and ion densities or the volume (or both) 
increase at these phase intervals. If these regions are associ- 
ated with the interaction sites between the compact object 
and the stellar wind this might imply that the stellar wind 
is restricted to these phases, e.g. it is a disk-like structure 



around the WR star as proposed in Fender et al. ( 1999 1 



For the COMPTT model the Comptonization normaliza- 
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Ep e [6.25,6.53] 
a Fe [0.14,0.38] 



bremsstrahlung, during a major radio flaring episode. From 



the definition of Koljonen et al. (2010 1 the change from the 



BN [0.70,1.68] 
kT bb/comp [0.82,1.19] 
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Figure 13. Phase-circle diagrams displaying the phase variation 
of fit parameters and bolometric X-ray flux (3.5-200 keV). Each 
of the four coloured diagrams present three fit parameters plot- 
ted on a circle, where the angle shows the orbital phase and the 
height of the curve shows the mean parameter value from stacked 
observations (the number of observations shown at the center di- 
agram for each phase bin) with the minimum value on the inner 
circle and the maximum value on the outer circle. In the brackets 
are the minimum and maximum values of a given parameter. The 
individual fits with parameters and their units can be found in 
Table[A] Upper left: absorption parameters: the covering fraction 
of the partial absorber (cov) and the optical depths of the ab- 
sorption edges of H- and He-like iron (tJ, and tJ, ). Upper right: 
iron line complex parameters: Equivalent width (EWp e ), cen- 
troid (Ep e ) and line width (ap e ). Lower left: Comptonization 
parameters: normalization (BN), soft seed photon temperature 
(kT bb i comp ) and the compactness of the radiation from soft seed 
photons (l s )- Lower right: bremsstrahlung parameters (normal- 
ization BrN and plasma temperature kT bb / br ) and bolometric 
X-ray flux (F^ i.) Center: the number of pointings contributing 
to a given phase bin. The concentric circles show the number of 
observations from (smallest circle) to 10 (largest circle). 



tion behaves exactly like the bremsstrahlung normalization. 
The only difference with the bremsstrahlung model comes 
from the inclusion of optical depth in COMPTT which pro- 
duces a maxima around phase 0.75. 

Based on the phase-folded principal components (see 
Fig. [6JI, we see that the first principal component has max- 
ima between phases 0.2-0.4 and 0.6-0.8, much in line with 
the belm normalization. However, a sharp peak is present 
at phase 0.75. A similar peak is observed also in the sec- 
ond principal component with additional maxima between 
phases 0.0-0.2 and around phase 0.5. 



5.3 Location of the Comptonizing region 

X-ray/radio state changes seem to be driven by the in- 
terplay of two spectral components, Comptonization and 



FHXR state to the FIM state signals softening of the X- 
ray spectra. This is coupled to a decrease of the radio flux 
density as the major radio flare fades (see Fig. [l] bottom 
left panel and Table [2j . The Compton upscattering seem 
to be more efficient during the FHXR state. Fig. [5] shows a 
stronger increase in the first principal component (which is 
coupled to the Comptonization normalization) during this 
state than for the second principal component (which is 
coupled to the bremsstrahlung normalization). Also, Fig. [7] 
shows that the X-ray flux attributed to Comptonization de- 
creases and the X-ray flux attributed to bremsstrahlung in- 
creases when the source moves from the FHXR state to the 
FIM state. Possible interpretations of this behaviour are: (1) 
the Comptonizing region is located in the jet and it is moving 
away from the pool of seed photons. Close to a major radio 
flare peak, most of the seed photons see the Comptonizing 
region and produce a harder spectra. When the Comptoniz- 
ing region moves away and/or dilutes as the major radio 
flare fades more seed photons arrive to us unscattered. (2) 
The increase of the seed photons (possibly corresponding to 
the increase of the stellar wind) cools the Comptonization 
region at the base of the jet so that there is a decrease in 
Comptonization and increase in unscattered seed photons. 
This would also serve as a mechanism to quench the jet ejec- 
tion episode. In essence these two scenarios can be dubbed 
as the distant location and the close location of the Comp- 
tonizing region, respectively. 

While it is possible that the high energy electron popu- 
lation originates close to the compact object it seems more 
probable that it originates further out, e.g. in the jet, be- 
cause of the low optical depth for the Comptonized emission. 
When using the above best-fitting model, we found that the 
Thomson optical depth for the Comptonization approaches 
its lowest possible value (in our model this corresponds to 
Tr=0.5). A close origin is also not favourable if the system 
inhabits a WR companion star, as its dense stellar wind most 
likely enshrouds the compact object causing extra scattering 
of the photons originating close to the compact object. 



5.4 Comparison to other microquasar systems 

In addition to Cyg X-3, similar thermal, hot components 
have been identified in other microquasar/XRB systems. 
Titarchuk fc Seifina| ( |2009| > found that eight IS/HSS spec- 
tra from GRS 1915+105 required a "high-temperature black 
body-like" profile; these were somewhat related to epochs of 



radio activity. Also, Mineo et al. (20121 found that spectra 
in the "heartbeat" state of GRS 1915+105 can be fitted with 
models including a black body component with colour tem- 



perature 3-6 keV. Seifina & Titarchuk (20101 found 24 IS 



spectra from SS 433 during radio outburst decay that had 
a strong black body-like component with colour tempera- 
ture 4-5 keV. In addition, the spectra from the black hole 
soft X-ray transients GS 2000+25, GS 1124-68 and XTE 
J1550— 564 were successfully fitted with an additional ther- 
mal Comptonization component with colour temperature 2- 
4 keV and high optical depth r ~ 5 ( |Zycki et al.|2001[ >. This 
prompts the question: is this thermal component something 
intrinsic to black hole systems/microquasars? Is the emis- 
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sion mechanism bremsstrahlung or thermal Comptonization 
as it appears to be in Cyg X-3 or something else? 

5.5 Origin of the thermal component 

In the papers mentioned above a multitude of different ori- 
gins were presented for the thermal component. Here we 
briefly review them and consider a few extra options as well. 
Since rigorous testing of the individual origin is beyond the 
scope of this paper, this discussion will be somewhat mini- 
mal. 

Gravitationally redshifted annihilation line? In 



Titarchuk & Seifina ( 2009 1 the thermal component was in- 



terpreted as a highly gravitationally redshifted annihilation 
line. However, a clear view close to compact object would 
also result in a better view of the accretion disc which is not 
observed, based on our modelling, at the same time with the 
thermal component. Rather, it seems that the impact of the 
disc is more prominent when the thermal component is at 
minimum. 

Hot disc? One option is that the accretion disc is 
partly heated by the hard X-rays (thus reaching colour tem- 
peratures higher than the usual ~ 1 keV). Physically, this 
would represent e.g. a patchy corona above the disk with 
hard X-rays arising from magnetic reconnection. However, 
based on the fluxes this scenario is not favourable. The hard 
X-ray flux ranges between fio-200k c v ~ 0.7 — 6.3 (in units 
of 10 -9 erg cm -2 s _1 ) while the thermal component flux 
(as modelled by thermal bremsstrahlung) ranges between 
F th ~ 1.0 - 7.9 with F th > Jio-aookeV in nearly all RXTE 
pointings. Also, we find the reflection component to be very 
small or non-existent and thus, reprocessed disc emission is 
not likely the origin of the thermal component. 

Three-phased accretion? This scenario is based on 



the idea ( |Zycki et al. 2001 1 that a transitionary warm plasma 
region could form between the cold accretion disc and the 
hot inner flow/corona. This would then supply seed photons 
for Comptonization in the corona. However, it is hard to 
envision how the orbital changes are incorporated into this 
scenario. 

Thermal jet? In a similar case to the previous point, 



Memola et al. (2002) studied a thermal plasma component 



in the jet. According to this study the innermost jet can 
reach temperatures above 1 keV up to 1 MeV, thus suggest- 
ing that thermal X-ray emission might originate from this 
region. The observations of the thermal components during 
intermediate states, during which strong jets are observed 
(especially for Cyg X-3), supports this scenario. However, 
the evolution of the principal components shows that while 
the major radio flare decays the thermal component becomes 
more prominent, which would point to a non-jet origin. 

Hybrid electron distribution? By using belm we 
assume a hybrid electron distribution in the coronal elec- 
trons ab initio motivated by the successful use of hybrid 
models in fitting Cyg X-3 spectra (e.g. |Koljonen et al.|2010| 
for COMPPS). Interestingly, COMPPS was found to success- 
fully fit spectra of GRS 1915+105 ( |Mineo et al.|2012[) XTE 
J 1550- 564 and some from GS 1124-68 ( [Zycki et al.|20"oTj ). 
However, as portrayed by the PCA results and the unsuc- 
cessful spectral fits using just belm, at least for Cyg X-3, 
the need for additional thermal component is evident. 



the presence of a Thomson-thick, low-temperature plasma 
cloud surrounding the compact object in Cyg X-3 in order 
to explain the lack of high frequencies in the power spec- 
tra and the peculiar hard state X-ray spectra with ~ 30 
keV cut-off by Compton downscattering. This cloud could 
be the result of the collision between the WR stellar wind 
and the compact object. The plasma parameters found in 
Zdziarski et all (120091) were kT ~ 3 keV and r ~ 7. For 



these parameters the thermal bremsstrahlung emission be- 
comes a substantial source for photons which get amplified 
by Comptonization in the plasma cloud. The orbital changes 
in the thermal component could arise if the wind is asym- 
metric and restricted to certain phases along the orbit. 



5.6 Synopsis 

Based on the above discussion we propose a scenario which 



explains the observed behaviour of the system (see Fig. 14 1 



Scattering cloud? Zdziarski et al. (20091 argued for 



(i) The principal emission process responsible for most 
of the variability in the X-rays during jet ejection events 
is the Compton scattering of soft seed photons from two 
distinct thermal populations by the high-energy electrons 
confined most likely in shocks in the relativistic jet. 

(ii) One source of these thermal populations is a ther- 
mal plasma cloud around the compact object, but which 
is restricted to phases 0.0-0.25 and 0.5-0.75. We interpret 
this (inspired by Zdziarski et al.|2009 \ as an interaction be- 
tween the compact object and the stellar wind. Therefore 
the stellar wind could be formed around the WR star as a 
disc. This is in line with the studies that show that approx- 
imately 15% of WR stars have anisotropic winds, which are 
most likely caused by equatorial density enhancements pro- 
duced by high rotation rates ( |Harries et al. 19981. In this 
plasma cloud the bremsstrahlung emission process becomes 
important and provides soft seed photon population that 
Compton upscatter in the jet. 

(iii) Outside the disc-like wind (provided that the com- 
pact object does not orbit the WR star exactly in the plane 
of the wind) this interaction is much weaker. The disc-like 
wind provides higher absorption when the compact object 
is under the disc along the line of sight (phase 0.25-0.5) 
and lower absorption when above the disc (phase 0.75-0.0). 
However, this is more a matter of going from a high density 
to a lower density part of the wind. But even in the low 
density part of the wind the density still needs to be rather 
high in order to wipe out the timing signature above 0.1 Hz 
( Koljonen et al.|2011| |Axelsson et al.|2009 >. 

(iv) During times outside the disc-like wind, a second 
source of soft seed photons arises from the accretion disc. 
Hints for glimpsing the underlying disc is obtained by the 
increasing X-ray variability during these phases. [Koljonen et| 
ak] ( |20li] ) found that during phases ~0.3 and ~0.7 the X-ray 
variability peaked above what could be expected based on 
the orbital modulation. Also, on two occasions a mHz QPO 
was detected around phase ~0.3. In addition, based on the 
spectral fits on the RXTE data the phase interval 0.3-0.4 
shows a minimum in the soft seed photon temperature as 
well as the centroid of the gaussian line complex approaching 
6.4 keV. 
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Figure 14. a) A sketch depicting the geometry of the system, with the Wolf-Rayet companion surrounded by a disc-like stellar wind and 
the companion object orbiting it (shown in four different orbital phases). Depending on the orbital phase the compact object is either 
inside or outside the stellar wind, b) The two different soft seed populations arise from the accretion disc and c) from the interaction of 
the compact object with the stellar wind. 



6 CONCLUSIONS 

In this paper we have studied in detail the X-ray behaviour 
of Cyg X-3 during a major radio flaring event. We have 
shown that additional information is obtained by using Prin- 
cipal Component Analysis on fast X-ray spectral data (~ 1 
minute), which helps to discriminate between different spec- 
tral models fitted to the time-averaged spectra (~ 3 ks). We 
have found that most of the variability is caused by two 
principal components whose evolution based on spectral fits 
is best reproduced by a Comptonization component and a 
bremsstrahlung component that most likely arises from the 
interaction between the compact object and the stellar wind 
from the Wolf-Rayet companion. Following the fit parameter 
evolution over the orbit, it appears that the bremsstrahlung 
component is restricted to certain phases that are opposite 
to each other, which could indicate that the stellar wind 
is formed as a disc-like wind around the companion star. 
Therefore, during orbital phases when the compact object 
is outside the wind we get a glimpse deeper into the directly 
of the accretion disc. This is also supported by the findings 
of increased variability and QPOs in the X-ray lightcurves 
by Koljonen et al. (20111. Thus, we propose a scenario for 



Cyg X-3 where the soft seed photon population arises from 
two different origins: from bremsstrahlung and accretion disc 
that in turn get Comptonized further down in the jet. 
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18 K. I. I. Koljonen et al. 

Table Al. phabs X pcfabs X edge(l) X edge(2) X (gauss + belm + bremss) 



Component 


Parameter 




Unit 


N2. 1 


2 


N2- 3 


N° 4 


N2. 5 


N2. 6 


phabs 


nH 




10 22 


3/ 


3/ 


3/ 


3/ 


3/ 


3/ 


pcfabs 


nHi 




10 22 


12/ 


12/ 


12/ 


12/ 


12/ 


12/ 




cov 






n O-1+0.04 

U-»l_o. 04 


0-86lH 3 4 


0.89l°;°04 




0.89ir 04 


64+ 04 
U - D *_0.04 


edge(l) 


E S?e n 




keV 


9.28/ 


9.28/ 


9.28/ 


9.28/ 


9.28/ 


9.28/ 




^-(i)i? e n 






17+ - 04 

U-J -' -0.04 


13+ 004 
u — 0.04 


n 17+0-03 
u ' x ' -0.04 


10+0.03 

u - ±0 — 0.04 


n 17 +0.03 

u ' x ' —0.04 


12+ 004 
u ' lz -0.04 


edge(2) 


TTiion 




keV 


8.83/ 


8.83/ 


8.83/ 


8.83/ 


8.83/ 


8.83/ 




^(2)i? c n 






14+ - 04 


12+ - 04 

U - iZ -0.04 


10+ 004 

u ' lu -0.03 


14+ 004 
"•^-o.os 


10+ 004 

u - lu -0.03 




gauss [Fe] 


EW Pe 




eV 


393 +2 4 


OS r; + 33 

2»a_ 30 


344+ 30 


382 + f 4 


299l| 


367lf 9 




Epo 




keV 


6 49+ - 02 


6 47+0.03 
—0.03 


6 43+0-03 
u " io -0.03 


6 44+ - 02 
"■^ -0.02 


6 43+0 03 
u -^ J -0.03 


f, 4 o+0.02 

b - 4S -0.02 




"Fe 




keV 




n 9s +0.08 
U - zS -0.09 


0.29+_° f 6 


24+ 05 
U - Z4 -0.06 


29+ - 06 
U - ZM -0.06 


14+ - 06 
u ' 1 ^-0.14 


BELM 


norm 








1 ,fi+ - 06 


l.5ll° ; 2 2 


i iq+0.12 
±,±o _0.13 


-i +0.20 
±mO1 —0.20 


1 oi+O.ll 
± - o± _0.11 




T 






0.5/ 


0.5/ 


0.5/ 


0.5/ 


0.5/ 


0.5/ 




1b 






100/ 


100/ 


100/ 


100/ 


100/ 


100/ 




Is 






i3.9tt;Si 


21.43 + 2'5 7 


20.92+ 2 -® 2 


20.531JS 


20.92+ 2 '® 2 


17.99+ 2 " 




kT bb 






i.05+;: 

— 0.07 


0.99+°°« 

—0.05 


0.82+°'°* 

— 0.07 


0.94+!!f 3 

—0.04 


0.81+°-°f> 

—0.07 


1.04+00 3 

— 0.04 


bremss 


norm 






z " y4l _0 68 


5.03±»-|? 


i oq+0.08 
i,di, _0 09 


1 12 +ai 


1 oq+0.08 
1 " J ^_0 09 


1.56±°-» 




kT bb 




keV 


1.21+"'' 2 

—0.14 


4.15+° : l 8 ) 

—0.12 


9.47+° : f, 

— 0.55 


9.09+° :7 i 

— 0.64 


9.47+?, : ^ 

—0.55 


7.34+0'^ 

— 0.59 


Confidence 


xL/ d -°- f - 






0.81/57 


0.84/63 


0.97/65 


1.06/65 


0.63/68 


0.86/63 


Fluxes 


pa 

bol,abs 


io- 9 


erg cm -2 s — 1 


7.4 


13.1 


10.8 


8.9 


14.6 


11.4 




FBELM/Fbol 






0.56 


0.60 


0.58 


0.60 


0.59 


0.49 










0.39 


0.37 


0.39 


0.36 


0.37 


0.48 




FLine/Fbol 






0.04 


0.04 


0.03 


0.04 


0.03 


0.04 


Component 


Parameter 




Unit 


N° 7 


N2. 8 


N2. 9 


N2. 10 


N2. 11 


N2. 12 


phabs 


nH 




10 22 


3/ 


3/ 


3/ 


3/ 


3/ 


3/ 


pcfabs 


nHi 




10 22 


12/ 


12/ 


12/ 


12/ 


12/ 


12/ 




cov 






90+ - 05 
u - au -0.05 


n 71+0.05 
u - ' i -0.05 


Q2+ - 03 

u - t ' -0.03 


g+0-02 


Q2+ - 03 
-0.06 


n 01 +0.06 
U ' 8i -0.03 


edge(l) 


E S?e n 




keV 


9.28/ 


9.28/ 


9.28/ 


9.28/ 


9.28/ 


9.28/ 










17+ - 04 
u ' -0.03 


14+ - 04 
U -- L *_0.04 


13+ 004 
u - lo _0.05 


21+ 02 
u - zl -0.02 


24+ - 02 
u - z ^_ 0.03 


n 90+0.02 

u - zo _0.02 


edge (2) 


E i? c n 




keV 


8.83/ 


8.83/ 


8.83/ 


8.83/ 


8.83/ 


8.83/ 










n+ 03 

u - lx -0.03 


13+ 004 


n 13+0-05 
U ' i,3 -0.05 


o.iitS:E 


12+ 003 

U - iZ -0.03 


10+ 003 

u ' lu -0.05 


gauss [Fe] 


EW Fe 




eV 


370+ 22 


302l 22 


422+ 47 
^ zz -39 


3101*? 


269+« 


389+3° 




E Fc 




keV 


6 42+ - 03 
u,y: -0.03 


fi ,,0+0.02 
b - 4S -0.02 


6 41+ 03 


6 42+ - 03 
°- £tz _0.04 


6 43+ - 04 


6 47+0-03 

»-1'_0.03 




fp e 




keV 


30+ 06 


14+ - 07 


4+ 06 
u - -0.06 


34+0.08 
u -°^-0.08 


0.34+00 9 


29+ 08 
u - za -0.08 


BELM 


norm 






1 26+ 019 


i 07 +0.05 
- L - U ' —0.04 


i 13+O.H 
— 0.08 


1 09+ - 05 
—0.07 


0.981^1 


0-87±°f 6 




r 






0.5/ 


0.5/ 


0.5/ 


0.5/ 


0.5/ 


0.5/ 




1b 






38.59/ 


100/ 


100/ 


100/ 


100/ 


100/ 




l s 






25.6915-62 


99+4.63 
zz -5.06 


12.64t 3 ;^ 


10.821 2 ; 93 


15.18_ 4 ; 25 


1 1 Q4+ 4 -0i 

ll.M1_ 1 g4 




kT bb 






n 77 +0.06 
u -' '-0.09 




95+ - 07 


I 14+0.09 


> 1.13 


> 1.15 


bremss 


norm 








3.19+ii 9 


4 43+0.76 
^•^°-0.74 


7 04+0.43 
' ,z -0.49 


19 9 +1.94 
lz - -4.69 


5 gg+18.66 
°- DM -0.64 




kT bb 




keV 


q 71 +0.39 
a - ' -0.38 


4 41+0-49 


o 9 o+0.08 
°- zo -0.09 


2 91+ 08 
Z - M1 -0.14 


9 1+0.58 
z -0.14 


2 61+ 24 
z - D1 -0.49 


Confidence 


X r 2 ed/ d -°- f - 






1.37/67 


0.92/58 


0.82/59 


1.10/61 


0.89/60 


0.78/53 


Fluxes 


pa 

bol,abs 


1Q -9 


erg cm -2 s — 1 


10.7 


10.0 


7.4 


8.1 


8.1 


6.4 




FBELM/Fbol 






0.62 


0.61 


0.58 


0.53 


0.63 


0.60 




FBremss/Fbol 






0.34 


0.36 


0.37 


0.43 


0.34 


0.35 




FLine/Fbol 






0.04 


0.04 


0.05 


0.04 


0.04 


0.05 



Frozen in the fits. 

a Absorbed, bolometric flux (3.5-300 keV) of the model normalized to the PCA data 
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Table Al - continued 



Component 


Parameter 




Unit 


N- 13 


N- 14 


N- 15 


N- 16 


N- 17 


N- 18 


phabs 


nH 




10 22 


3/ 


3/ 


3/ 


3/ 


3/ 


3/ 


pcf abs 


nHi 




10 22 


12/ 


9/ 


9/ 


9/ 


9/ 


9/ 




COY 






QO+ - 03 
u ' az -0.04 


> 0.90 


0.90tg$ 


n 47+0-07 
u -^< -0.10 


0-95^. 3 


0.95«.oi 


edge(l) 


E S?e n 




keV 


9.28/ 


9.28/ 


9.28/ 


9.28/ 


9.28/ 


9.28/ 




^(l)i?e n 






n 9 o+0.02 
"■^-0.02 


09+ 04 

u - ua -0.04 


14+ 004 

U ' i ^-0.05 


oo+ - 01 


29+ 02 
u - za -0.02 


n 10+0. 03 
U - iS -0.03 


edge(2) 


E i? c n 




keV 


8.83/ 


8.83/ 


8.83/ 


8.83/ 


8.83/ 


8.83/ 




^(2)i? c n 






n in+ 03 

u ' lu -0.03 


17 +0.05 
U - i ' -0.05 


o.09±g;» 


n 2^+0.02 

u - zo -0.02 


n+0.01 
u -o 


o,4+ - 04 

u - u4 -0.03 


gauss [Fe] 


EW Pe 




eV 


341 + 34 




305± 3 


546l 23 


416l 2 2 9 6 


ozo -105 




E Pe 




keV 


R ak+0.03 
D -*°_0.03 


fi 42+0- 03 
D -^ z -0.03 


a cc+0.03 
D - oo -0.03 


6.26±° : {g 


6-63l°'°3 


R Ql +0.07 
o - J1 -0.06 




l^Fe 




keV 


94+0.07 
u ' z4 -0.08 


21+ 009 


0.26+. ,' 08 


1 ,+°- 06 


40+ 05 

u '^ u -0.05 


74+0.09 
u -'^-0.12 


BELM 


norm 






i on+0.06 
i.oa_Q 05 


n 70+O.04 
u - lo -0.04 


90+ 002 

u ' au -0.03 


n 24+0- 01 
u - Z4 -o.oi 


n rc+0.03 

u - DO -o.oi 


n q 2 +0- 03 
u ' az -0.O3 




T 






0.5/ 


0.5/ 


0.5/ 


0.5/ 


0.5/ 


0.5/ 




1b 






12.63+S; 92 


29.07+. 42 ,, 72 


49 15+ 49 - 34 


76 4+123.6 
'"•^-64. 88 


21.18+ 22 ; 7 


20 01 + 11 - 73 
/u.ui_ g 4g 




Is 






< 10.55 


I8.l8l 6 5 ;f 8 


22.71+gg 7 


10+ 4 ' 62 


lot 40 


10 +1.95 




kT bb 






1.30+°?! 

— 0.11 


1.03+S-S5 

— 0.04 


1.06+°/ n 3 

— 0.05 


i.i6 + °-°S 

— 0.05 


> 1.34 


1 10+0. 08 
— 0.07 


bremss 


norm 






7.58+1-™ 


gn+0.89 


1 cc + 1.31 

1-661 1 06 


79+0-05 
u - ' a _o 25 


6 57+ - 17 
D -°'_0 65 


K ,4+0-42 
°- o4 _0 42 




kT bb 




keV 


2.13+Of, 

— 0.15 


2 qo+0-20 
—0.13 


— 0.38 


ok+0.32 

— 0.15 


2 00+O.IO 
— 0.04 


9 qo+0.12 
— 0.15 


Confidence 


X 2 Icd /d.o.{. 






0.97/ 72 


0.96/58 


1.07/ 64 


1.22/49 


1.24/ 64 


1.49/ 69 


Fluxes 


pa 

bol , abs 


10 -9 


erg cm -2 s — 1 


7.9 


5.8 


7.1 


2.0 


5.5 


6.2 




FBELM/Fbol 






0.74 


0.64 


0.74 


0.53 


0.55 


0.65 




^Bremss /-^bol 






0.22 


0.33 


0.23 


0.41 


0.40 


0.30 




FLine/Fbol 






0.03 


0.03 


0.04 


0.06 


0.04 


0.06 


Component 


Parameter 




Unit 


N- 19 


N- 20 


N- 21 


N- 22 


N- 23 


N- 24 


phabs 


nH 




10 22 


3/ 


3/ 


3/ 


3/ 


3/ 


3/ 


pcf abs 


nHi 




10 22 


9/ 


9/ 


9/ 


9/ 


9/ 


9/ 




cov 






0.95+° 01 


qo+0-02 
u - ao -0.04 


0.95±° 01 


0.80+°- 95 5 


n yo+0.05 
u - ,o -0.06 


95+° 
u - Md -0.02 


edge(l) 


T?ion 




keV 


9.28/ 


9.28/ 


9.28/ 


9.28/ 


9.28/ 


9.28/ 




^(l)S?o n 






15+0- 02 
u - io -o.oi 


20+ - 04 
U - ZU -0.04 


11+ 004 


02+ - 04 
u - uz _0.02 


14+ 004 

u - 14 -0.04 




edge (2) 


E i? c n 




keV 


8.83/ 


8.83/ 


8.83/ 


8.83/ 


8.83/ 


8.83/ 




^(2)i? e n 






u ' lo -0.02 


O.Q8+.°;°05 


n 14+005 
u - 14 -0.05 


n OO+0.03 
u - zz -0.04 


o-i5l :°ol 


n ia+0.04 


gauss [Fe] 


EW Pe 




eV 


240+23 
Z4J _20 


4301^ 


364+ 42 

—35 


343+ 2 » 


486+^ 


323+ 39 




Ef 6 




keV 


f? 90+0.03 
"■ zo -0.04 


R cq+0.02 


R 44+0.03 
D - 44 -0.03 




R 4O+0 2 
°' 49 -0.01 


6-40l°;° 3 








keV 


n 2c:+0.07 
u ' zo -0.09 


n 99+O.O6 
u - zz -0.08 


30+ 07 
u - JU -0.08 




02+ 08 
u - uz -0.02 


OO+0.07 


BELM 


norm 






1 02+ 03 
i - UZ -0.04 


c-y+0.05 
U -°'-0.05 


60+ 03 

u ' ou -0.04 


0.8<« 


n 09+O.02 
u - oz -0.02 


-. 01 +0.05 
1 - ol -0.06 




T 






0.5/ 


0.5/ 


0.5/ 


0.5/ 


0.5/ 


0.5/ 




1b 






9 o 4+14.86 
zo - -8.26 


100/ 


100/ 


100/ 


100/ 


100/ 




Is 








16 21+ 319 


20.59+2 24 


99 c;^+ 2 - 8 
zz.oo_ 2 6 3 


10 47+ 3 - 14 
1Ui4 '-0.47 


21.74+2 7 9 




kT bb 






1 19+ 012 
0.07 


^ + -° °ot 


1 01+ 003 
J - ul -o.o3 


1 01+ 004 
1 - ul -0.04 


1 1 7 +0.06 
-0.06 


OO+0.03 

u - ay -0.03 


bremss 


norm 






10 99+ ' 93 


0.86+?;?° 


75+0- 15 
U -'°-0.12 


1 64+ 033 


0.87t° .f 9 


1 ss+0- 51 

l-»»_0.42 




kT bb 




keV 


41 +0.06 


6 45+ - 76 
D-^O-0.49 


R ,0+0 63 
o - oa -0.46 


5 4^+0.47 
O - 4O -0.34 




r ci +0.60 
°- o± -0.40 


Confidence 








1.29/75 


1.23/57 


1.13/60 


0.84/61 


0.91/56 


1.01/62 


Fluxes 


pa 

bol, abs 


lO" 9 


erg cm -2 s — 1 


8.4 


4.6 


4.7 


7.3 


2.8 


10.1 




FBELM/Fbol 






0.52 


0.60 


0.66 


0.60 


0.48 


0.67 




FBremss/Fbol 






0.45 


0.36 


0.30 


0.35 


0.47 


0.29 




FLinc/Fbol 






0.03 


0.04 


0.04 


0.04 


0.05 


0.04 



Frozen in the fits. 

a Absorbed, bolometric flux (3.5-300 keV) of the model normalized to the PCA data 
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Table Al - continued 



Component 


Parameter 


Unit 


25 


N2- 26 


27 


phabs 


nH 


10 22 


3/ 


3/ 


3/ 


pcf abs 


nHi 


10 22 


9/ 


9/ 


9/ 




cov 




n S5+ ' 05 

U -°°-0.05 


70+ - 06 

u -' u -0.06 


n s1 +0.05 


edge(l) 


pion 
Fe 


keV 


9.28/ 


9.28/ 


9.28/ 




T ( 1 )S?e 11 




o-io±8:St 


o-i6±S:g| 


0-17i^4 


edge(2) 


E i? c n 


keV 


8.83/ 


8.83/ 


8.83/ 




T(2)S? o n 




17+ - 04 

u - i ' -0.04 


13+ 005 




gauss [Fe] 


EW Fe 


eV 


316+JI 




299+29 




Ep c 


keV 


6 47+ - 02 

"■^'-0.03 


6 51+ 02 


6 47+0.03 

"•^'-0.03 




<7p e 


keV 


20+ ' 06 
u - zu _0 08 


0.13±g-°| 


31+ 07 
U - J± _0 07 


BELM 


norm 




-, --+0.12 
L - 1 °_o 11 




-, 77 +0.16 
A - ' ' _o 16 




r 




0.5/ 


0.5/ 


0.5/ 




Id 




38.59/ 


100/ 


100/ 




Is 




18.07+ 2 ; 2 , 


-,0+7.92 
18 -4.41 


17.82+ 2 ;« 




kT bb 




1 n4 +0 ' 04 


-, 19 +0.05 
1 - lz -0.05 


o.9»±8:SS 


bremss 


norm 




o nl +0.53 


o.90tg;i; 


1 99+ 034 

±.33_g 27 




kT bb 


keV 


6 27+ ' 46 
D,Z '-0.39 


9.65+jj 3 


7 n-,+0.72 


Confidence 


X? ed /d.o.f. 




1.27/64 


1.02/53 


0.89/64 


Fluxes 


pa 

bol.abs 


10 — 9 erg cm -2 s — 1 


15.0 


6.7 


13.6 




FBELM/Fbol 




0.58 


0.52 


0.63 




^Bremss /Fbol 




0.38 


0.44 


0.34 




FLine/Fbol 




0.03 


0.05 


0.03 



f Frozen in the fits. 

a Absorbed, bolometric flux (3.5-300 keV) of the model normalized to the PCA data 



pcfabs X edg 


e(l) X edge(2) X 


(linemodel + belm 4 


- bremss) 






Component 


Parameter 


Unit 


N° 1 


N° 2 


N° 3 


phabs 


nH 


10 22 


o 99 + 0.12 
^•^-O.ll 


4 18+ ' 09 

* ,lo -0.09 


3.ooiS:" 


pcfabs 


nHi 


1Q 22 


in 44+ - 52 

iu.^_ 4g 


16 2+ 08 


8.98±° £ 




cov 




n 70+0 01 

u - ,M -0.01 


n 70+0.01 
u - '^-O.Ol 


0.85ir i 


edge(l) 


Ei?" 


keV 


9.28/ 


9.28/ 


9.28/ 








n 15+ 002 

u ' io -0.03 


14+0-03 


09+ - 03 
u - uy -0.03 


edge (2) 




keV 


8.83/ 


8.83/ 


8.83/ 




^(2)l?e n 




n+ 03 


10+ 002 

u,lu -0.02 


1T+0.03 
V- 11 -0.03 


BELM 


norm 




1 02+ 03 
1 - u -0.03 


1 39+ 004 


77+0-01 
u -' '-0.01 




r 




0.5/ 


0.5/ 


0.5/ 




1b 




> 122 


> 185 


loilg 




Is 




18 9+ 14 

1 00+ 02 
1 - uu -0.02 


23 fi+ ' 8 
0.98lg : S 


< 11.0 




kT bb 




-1 17 +0.04 
1-1 ' -0.03 


bremss 


norm 




-, 99 +0.09 
1,z -0.09 


i-02l° :S 


4 48+ - 14 




kT bb 


keV 




7 fic -+0.33 
'• D0 -0.31 


2 89+ - 04 
z -°»_0.05 


Confidence 


X?ed/ d -°- f - 




1.48/225 


1.57/ 291 


1.61/217 


Fluxes 


pa 

bol,abs 


10 — 9 erg cm" 2 s" 1 


10.9 


10.5 


5.7 




FBELM/Fbol 




0.62 


0.72 


0.55 




FBremss /Fbol 




0.34 


0.25 


0.41 




FLine/F bo l 




0.04 


0.03 


0.04 



f Frozen in the fits. 

a Absorbed, bolometric flux (3.5-300 keV) of the model normalized to the PCA data 
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